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The quantum phase transition in clean itinerant ferromagnets is analyzed. It is shown that soft 
particle-hole modes invalidate Hertz's mean-field theory for d < 3. A renormalized mean-field 
theory predicts a fluctuation-induced first order transition for I < d < 3, whose stability is analyzed 
by renormalization group techniques. Depending on microscopic parameter values, the first order 
transition can be stable, or be pre-empted by a fluctuation-induced second order transition. The 
critical behavior at the latter is determined. The results are in agreement with recent experiments. 



PACS numbers: 75.20.En; 75.45. +j; 64.60. Kw 

One of the most common and basic examples of a phase 
transition is the paramagnet-to-ferromagnet transition in 
metals, e.g., iron and nickel. These elements have high 
Curie temperatures, on the order of 1,000 K. Examples 
of itinerant ferromagnets with much lower Curie tem- 
peratures, on the order of tens of K, include MnSi [1], 
ZrZn2 [1,2], and UGe2 [3]. In the latter materials, the 
Curie temperature can be suppressed to zero by the ap- 
plication of hydrostatic pressure. This allows for an ex- 
perimental investigation of the ferromagnetic quantum 
phase transition, which takes place at zero temperature 
as a function of a non-thermal control parameter, in this 
case, pressure. Another example is Ni K Pdi_a; [4], where 
the control parameter is the nickel concentration. 

The ferromagnetic transition in metals was also the 
subject of the earliest theoretical studies of quantum 
phase transitions. In an important paper, Hertz [5] con- 
cluded that the critical behavior should be mean field-like 
in all dimensions d > 1. This is because, in quantum sta- 
tistical mechanics, statics and dynamics are coupled. As 
a result, the quantum phase transition in d dimensions 
is related to its classical counterpart in d + z dimensions, 
with z the dynamical critical exponent. Since simple the- 
ories suggest z = 3 for clean itinerant ferromagnets, and 
since the classical Heisenberg ferromagnet has an upper 
critical dimension <i+ = 4, this argument suggests d£ = 1 
for the quantum transition. 

It is now known, mostly through studies of the cor- 
responding problem in the presence of quenched disor- 
der, that the above argument is in general not correct. 
The basic physical reason is the existence of soft modes, 
particle-hole excitations in the case of itinerant ferromag- 
nets, that couple to the order parameter and preclude the 
construction of a Hertz-type Landau-Ginzburg- Wilson 
(LGW) theory entirely in terms of the order parameter. 
These soft modes lead to time scales in addition to the 
critical one, and hence to multiple exponents z. This in 
turn leads to an instability of the mean-field fixed point 
that is not apparent in a power-counting analysis of the 



LGW theory. In the disordered case, the net result is 
an upper critical dimension d+ = 4 (instead of d+ = 
in Hertz theory), with non-mean field (and non-power 
law) critical behavior in the physically most interesting 
dimension d = 3 [6]. In the clean case, an analogous 
analysis of the stability of Hertz's fixed point shows that 
the upper critical dimension is rf+ = 3. In d = 3 a gen- 
eralized Landau theory predicts a first order transition 
due to an m 4 In m term in the Landau free energy, with 
m the magnetization [7]. For 1 < d < 3 this theory also 
predicts a first order transition. Sufficiently high temper- 
ature or disorder lead to an analytic Landau free energy 
and render the transition second order. 

This theoretical situation is at best in partial agree- 
ment with experiments. The existing theory predicts 
that the transition should always be of first order in 
sufficiently clean materials at sufficiently low tempera- 
tures [8]. If temperature or disorder drives it second or- 
der, then the predicted critical exponents are the classi- 
cal Heisenberg exponents, or the strongly non-mean field 
exponents of Ref. [6], respectively. Experimentally, the 
transition in MnSi and UGe2 at low temperatures is in- 
deed observed to be of first order [1,3], but in ZrZn2 it 
is of second order even in very clean samples at very low 
temperatures [2], and the same is true in NiPd [4]. Fur- 
thermore, the critical behavior in NiPd was found to be 
mean field-like to within the experimental accuracy, in 
good agreement with the predictions of Ref. [5] . This is 
suprising, given the above conclusion that Hertz theory 
cannot be correct in d = 3. 

In this Letter we show that the nature of the clean fer- 
romagnetic quantum phase transition is determined by 
physical effects that had not previously been recognized, 
and that taking these effects into account removes the 
above discrepancies between theory and experiment. We 
will use renormalization group (RG) techniques to show 
the following. (1) The first order transition can be un- 
derstood as a fluctuation-induced first order transition. 
(2) For certain microscopic parameter values the first or- 
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der transition is unstable with respect to a fluctuation- 
induced second order transition. (3) The upper critical 
dimension is rf+ = 3. (4) In the second order case, the 
critical behavior in d — 3 is given by mean-field expo- 
nents with logarithmic corrections, and in d < 3 it can 
be controlled by means of a 3 — e expansion. We will 
present our results first, and then sketch their derivation 
and explain their physical origin. 

If the bare value of the quartic coefficient in the Landau 
free energy is sufficiently small (in a sense to be specified 
below), one finds the first transition discussed in Ref. 
[7]. However, if the bare quartic coefficient is sufficiently 
large, the transition is of second order. In d = 3, the 
critical behavior is mean field-like with logarithmic cor- 
rections to scaling. Specifically, the paramagnon prop- 
agator in the critical regime in the paramagnetic phase 
has the form 



M(k, ft„) = l/(i + a(k) k 2 + \a n \/k) 



(la) 



where t is the dimensionless distance from criticality at 
zero temperature (T = 0), k is the wavenumber, and 
fi„ = 2irTn is a bosonic Matsubara frequency, k and 0„ 
have been made dimensionless by means of suitable mi- 
croscopic scales. The leading behavior of the coefficient 
a for small k is 



a(fc^0) oc (lnl/fc)- 1 / 26 



(lb) 



Such logarithmic corrections to power-law scaling can be 
conveniently expressed in terms of scale dependent crit- 
ical exponents. For instance, with b ~ 1/k a RG length 
scale factor [9], we can write a(k) k 2 oc k 2 ~ n , with a scale 
dependent critical exponent r\ given by 



rj = In In bl In b . 

1 26 1 



(2a) 



The correlation length exponent v, the susceptibility ex- 
ponent 7, and the dynamical exponent z can be directly 
read off Eqs. (1). The order parameter exponents (3 and 
S can be obtained from scaling arguments for the free 
energy. We find 



u = 1/(2 -rj) , z = 3-r 1 
7-1 , 0=1/2 , <5 = 3 



(2b) 



These exponents are defined as usual, i.e., £ <~ t v , O ~ 
T ~ £~ z , M ~ i" 7 , m ~ f 3 , and to ~ h 1 / 5 , with £ the 
correlation length and h an external magnetic field. The 
result for rj is valid to leading logarithmic accuracy; the 
values of 7, /3, and S, as well as the relations between rj 
and v and z, respectively, are exact. Finally, we define 
a specific heat exponent a by Cy oc T~ a at criticality. 
(This is a generalization of the usual definition of a at 
thermal phase transitions.) We obtain the exact relation 



a = — 1 — rf/z — In In 6/ In 6 



In d = 3 — e, the second order transition can be treated 
in an e-expansion. We find 



V 



-e/26 



a 



-l + (e-r])/z 



(3) 



The value for 77 is valid to one-loop order, the second 
equality is exact. The other exponents are still given by 
the exact Eqs. (2b) above. 

These results predict that the transition in 3 — d, to the 
extent that it is of second order, is characterized by mean- 
field exponents with logarithmic corrections. Within the 
accuracy of existing experiments, this is indistinguishable 
from Hertz's critical behavior. In addition to explaining 
why the transition is of first order in some materials, and 
of second order in others, our theory therefore provides 
an explanation for the fact that the observed critical be- 
havior in the second order case is mean field-like. We 
now sketch the derivation of the above results. 

The lesson learned from the disordered quantum fer- 
romagnetic problem [6] is the following: For a reliable 
analysis of the critical behavior it does not suffice to con- 
struct a LGW theory. Rather, in addition to the order 
parameter fluctuations, all other soft modes that couple 
to the latter must be kept explicitly and on equal foot- 
ing. Accordingly, the effective action should consist of a 
part depending on the order parameter field M, a part 
depending on the soft fermionic two-particle modes de- 
scribed by a field q, and a coupling between the two, 



A[M, q] = A M +A g + A 



M,q 



(4) 



Am is a static LGW functional (the dynamics will be 
provided by AM, q ), 



A M = j dx M.(x) [t - aV 2 } M(x) 



dxM 4 (x). (5) 



Here x = (x, r) comprises the real space position x 
and the imaginary time r. J dx = J dx. dr with 
(3 = 1/ksT. a and u are constants. 

The soft fermion field q originates from the composite 
fermion variables [10] 



Q12 = 



( -ipifihi -ipittn -V'itV^i 
-_V> u ;02t -i>ni>n -V'uV^i 

V -V'lT^T -V'lT^j. -V'ltV^l 



V>it^2T \ 
V>i|^2T ) 



(6a) 



(2c) 



Here the ip and ■tjj are the Grassmann-valucd fields that 
provide the basic description of the electrons, and all 
fields are understood to be taken at position x. The 
indices 1, 2, etc. denote the dependence of the fields on 
fermionic Matsubara frequencies to m — 2^T{ri\ + 1/2), 
etc., and the arrows denote the spin projection. A con- 
venient basis in the space of 4 x 4 matrices is given by 
r r ® Si (r, i = 0, 1, 2, 3), with r = s the 2x2 unit ma- 
trix, and 71,2,3 = — s i,2,3 = — ici,2,3, with the <7j the Pauli 



2 



matrices. The matrix elements of Q are bilinear in the 
fermion fields, so Q-Q correlation functions describe two- 
fermion excitations. In a Fermi liquid, the Q fluctuations 
are massive and soft, respectively, depending on whether 
the two frequencies carried by the Q field have the same 
sign, or opposite signs. We thus separate the Q fluctua- 
tions into massless modes, q, and massive modes, P, by 
splitting the matrix Q into blocks in frequency space [10], 

Qnm(x) = 9 (ran) P„ m (x) + 0(n)9(-m) <? nm (x) 

+e(-n)e(m)g+ m (x) . (6b) 

In what follows, we will incorporate the frequency con- 
straints expressed by the step functions into the fields P 
and q, respectively. That is, the frequency indices of q 
must always have opposite signs. 

The massive modes can be formally integrated out to 
obtain an effective action for the soft modes, q nm - The 
Gaussian part of the fermionic action has the form 

Aq = ~T~ f dxdy tr (9i2( x ) r i2!34( x -y)<?L(y)) i 

1,2,3,4 

(7a) 

where tr traces over the matrix degrees of freedom of 
Eq. (6a). As we see from Eq. (6a), the q propagator de- 
scribes particle-hole excitations, which in a clean electron 
system have a ballistic dispersion relation, i.e., the fre- 
quency scales linearly with the wavenumbcr. The vertex 
function r( 2 ) in momentum space therefore has the form 

(k + GH\n^ 2 \) . (7b) 

G and H are model dependent coefficients. A spin-singlet 
interaction amplitude can be included in the model, but 
will not be of qualitative importance for our purposes. 

The coupling term Am,q originates from the linear cou- 
pling between M and the electron spin density, which 
can be expressed in terms of Q by means of Eq. (6a). 
It is obvious that integrating out the massive field P 
will result in terms that couple M with all powers of 
q, Au, q = Am-q + Au-q 2 + • • •• Let us define a matrix 
magnetization field B(x) by 

Bia(x) = fa ® Si) (-r +1 *B(x) , (8a) 

i,r 

with components 

*B 12 (x) = $>„, ni -n 2 [M*(x) + (-) r+1 Mi n (x)] , (8b) 

n 

The first term in that series can then be written 

A M - q = ciT 1 / 2 J dx tr (B(x) g(x)) . (9a) 
The second one has the overall form 



A M - q i » ^T 1 ' 2 J dx tr (B(x) g(x) gt (x)) . (9b) 

Its detailed structure will be derived elsewhere [11]. c\ 
and C2 are model dependent coupling constants. 

The Gaussian field theory defined by the terms bilinear 
in M or B and q is easily diagonalized in terms of the 
paramagnon propagator 

and the fermion propagator 

V{k,9, n ) = l/(k + GH\a n \) , (10b) 

We now subject the entire action to a RG analysis [12]. 
We employ a differential momentum-shell RG and inte- 
grate over all frequencies. With b the RG length rescaling 
factor, we rescale wavenumbers and the fields via 

k bk , (11a) 
M„(x) - &( d - 2 +">/ 2 M rl (x) , (lib) 
q nm (x) - fo( d - 2 +«/V im (x) . (11c) 

Here M„(x) is the temporal Fourier transform of M(a;), 
and r\ and fj are exponents that characterize the spa- 
tial correlations of the order parameter and the fermion 
fields, respectively. The rescaling of imaginary time, fre- 
quency, or temperature is less straightforward. We need 
to acknowledge the fact that there are two different time 
scales in the problem, namely, one that is associated with 
the critical order parameter fluctuations, and one that is 
associated with the soft fermionic fluctuations. Accord- 
ingly, we must allow for two different dynamical expo- 
nents, z and z, and imaginary time and temperature may 
get rescaled via either one of two possibilities, 

r^b- z T , T-,b z T , (lid) 
T^b-~ Z T , T — > b z T , (lie) 

How these various exponents should be chosen is dis- 
cussed below. 

Within this framework, Hertz's theory corresponds to 
a fixed point where f\ = z = 1, which makes G and 
H marginal, and rj — 0, z — 3, which makes a and c\ 
marginal. Power counting then shows that C2 is irrelevant 
for d > 1 if the time scale is given by the exponent z, but 
marginal for d = 3 and relevant for d < 3 if it is given by 
z. It is easy to find explicit diagrams, starting at one-loop 
order, where the latter is the case [11]. This establishes 
that the upper critical dimension, above which Hertz's 
fixed point is stable, is c?+ = 3. 

To deal with the situation in d < 3 we go to one-loop 
order. Perturbation theory, combined with power count- 
ing, shows that in d = 3 there are no logarithmic correc- 
tions to ci, C2, G, and t. Motivated by the disordered 
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case [6] , we will be looking for a fixed point where G and 
ci are marginal, which implies 77 = 1 and z + 77 = 3. As 
mentioned above, C2 can have two different scale dimen- 
sions, depending on which time scale enters. We require 
C2 with the fermionic time scale, corresponding to z, to 
be marginal, which yields z + 77 = f . This leaves 77 as 
the only independent scale dimension. The critical time 
scale makes c-i irrelevant with scale dimension — 1. For 
the remaining coupling constants to one-loop order we 
find the flow equations 

da/d\nb = -rja - A a /H , (12a) 
du/d\nb = -(2 + 77)77- A u cj/H , (12b) 
dH/d\nb = r]H + A H /(a + t) . (12c) 

Here c 2 is the irrelevant version of C2, with flow equation 

dc 2 /dlnb = -c 2 . (12d) 

A u > is a model dependent positive coefficient. An 
explicit calculation yields for the other two coefficients, 
Ah = 27 A a = 3Gc 2 /ir 3 , with c 2 the marginal incarnation 
of this coupling constant. The ratio An/A a = 27 is what 
determines the critical exponents. 

Solving the flow equations, Eqs. (12), at t = shows 
that u becomes negative at a finite scale provided that its 
bare value satisfies < A u (c 2 °^) 2 /A a . This results in 
a fluctuation-driven first order phase transition. Indeed, 
the generalized mean-field theory for this transition maps 
onto the one for the superconducting transition at T > 0, 
which is the canonical example of a fluctuation-driven 
first order transition [13]. However, if the opposite in- 
equality holds, then u remains positive at all scales and 
the transition is continuous. Notice that within strict 
perturbation theory the second term on the right-hand 
side of Eq. (12b) is constant, so the transition is always 
of first order. This is the RG version of the theory given 
in Rcf. [7]. We note that for sufficiently large t, perturba- 
tion theory is valid. If the first order transition predicted 
by perturbation theory occurs at sufficiently large t, it 
therefore is unaffected by the fluctuation effects. In or- 
der for the generalized mean-field theory of Ref. [7] to 
be controlled, the magnetization discontinuity must in 
addition be small. 

As we see, fluctuations can qualitatively change the 
prediction of perturbation theory and drive the transition 
second order. The mechanism for this is very similar to 
the fluctuation-driven second order transition in classical 
Potts models [14]. The point is that the renormaliza- 
tion of a negative loop correction to the w-flow equation 
can make this term go to zero for large scales sufficiently 
fast to keep u positive, even if a scale independent loop 
correction would lead to a negative u. In the case of a 
second order transition, the critical exponents 77, v 1 and 
7 can be obtained by solving Eqs. (12) and substituting 
the result for a(ln b) in the paramagnon propagator, Eq. 



(10a). In d = 3 this procedure is tricky since it leads to 
scale dependent exponents; technical details will be re- 
ported elseswhere [11]. The results are given in Eqs. (2). 
In d = 3 — e the procedure is straightforward and leads 
to 77 as given in Eqs. (3), and to Eqs. (2b). a is most eas- 
ily obtained from the Gaussian free energy by replacing 
the Gaussian paramagnon propagator by the critical one. 
Differentiation with respect to T yields the specific heat. 
(3 and 6 are most easily obtained from scaling arguments 
[11], taking into account that u is a dangerous irrelevant 
operator for the magnetization [12]. 
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